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Abstract 

In this letter we present a new analysis of the torques acting on the accreting millisecond X-ray pulsars 
SAX J1808. 4-365 and XTE J1814-338, and show how our results can be used to constrain theoretical 
models of the spin evolution. In particular we find upper limits on any spin-up/down phase of XTE 
J1814-338 of |i>| < 1.5 X 10~^''Hzs~^ at 95% confidence level. We examine the possibility that a 
gravitational wave torque may be acting in these systems and suggest that a more likely scenario is 
that both systems are close to spin equilibrium, as set by the disc/magnetosphere interaction. 
Subject headings: pulsars: individual (XTE J1814-338, SAX J1808. 4-3658) — gravitational waves 



1. INTRODUCTION 

Accreting millisecond X-ray pulsars (AMXPs) are one 
of the best laboratories to test our understanding of 
strong gravity and of matter at extreme densities. In 
these systems a neutron star (NS) in a binary accretes gas 
from a less evolved low mass donor companion. During 
the accretion process angular momentum is transfered 
to the neutron star, which can be spun up to millisecond 
periods. 

This picture is not, however, without its complications. 
As more AMXPs and millisecond radio pulsars are dis- 
covered there appears to be strong evidence for a cutoff 
in the distribution of the spin rates at approximately 
730 Hz, which is well below the Keplerian breakup fre- 
quency for these objects (Chakrabarty et al. 2003; Pa- 
truno 2010). It thus appears that the torques acting on 
these systems are stopping them from spinning up to the 
breakup limit. A detailed analysis of the torques acting 
on an AMXP in a low mass X-ray binary (LMXBs) is a 
challenging problem, both from a theoretical and an ob- 
servational point of view. Observationally the problem is 
complicated by the presence of strong timing noise, that 
appears to be correlated with the X-ray flux in at least 6 
AMXPs (Patruno, Wijnands & van der Klis 2009). The- 
oretically, on the other hand, it is believed that the NS 
could also be spun down during accretion if the angular 
momentum transferred from the gas is not high enough, 
in what is known as the "propeller" phase (Illarionov & 
Sunyaev 1975). In fact recent simulations have shown 
that spin-down torques may be present for a wide range 
of accretion rates (Rappaport et al. 2004; D'Angelo & 
Spruit 2011). Furthermore the rapidly rotating NS may 
be emitting gravitational waves (GWs) that remove an- 
gular momentum from the system at a rate that is suffi- 
cient to halt the spin- up at a frequency of w 700 Hz (Bild- 
sten 1998). Finally, accretion torques that spin- up the 
NS act only during relatively short outbursts, whereas 
the NS is spun down by magneto-dipole torques during 
long quiescence periods. Therefore the AMXP spin fre- 
quency might reach a long-term equilibrium even if ac- 
cretion torques are present during outbursts ( Patruno 
2010). 

We present a new analysis of the torques acting on the 



AMXP XTE J1814-338 (referred to as XTE J1814) and 
use the results, together with those already obtained by 
Hartman et al. (2008, 2009) for SAX J1808.4-365 (SAX 
J1808), to constrain theoretical models of the spin evo- 
lution. 

2. OBSERVATIONS & DATA ANALYSIS 

2.1. SAX J1808.4-3658 

SAX J1808.4-3658 is a 401 Hz AMXP in a 2.01 hr bi- 
nary orbit (Wijnands & van der Klis 1998; Chakrabarty 
& Morgan 1998) that has undergone six outbursts since 
1996. Five outbursts have been monitored with the Rossi 
X-ray Timing F,xploi-er{RXTE). No significant spin- 
up/down episodes are detected during these outbursts, 
with upper limits of the order of Iz>| < 2.5 x 10"^'' Hz s~^ 
(Hartman et al. 2008, 2009). A long term spin-down is 
detected when comparing the constant spin frequencies 
measured in different outbursts. The measured long- 
term spin-down is i'sd — —5.5 x lO^^^Hzs^^, which 
has been interpreted as due to magneto-dipole torques 
induced by a NS magnetic field of '--^ 1.5 x 10®G. For this 
source we will use the results reported by Hartman et al. 
(2008, 2009) and will not perform any new data analysis. 

2.2. XTE J18U-338 

XTE J1814 was discovered in 2003 (Markwardt & 
Swank 2003), during an outburst that lasted nearly 2 
months. Only one outburst has been detected so far 
with an extensive coverage of RXTE whereas a tenta- 
tive recurrence time of 19 yr has been proposed based on 
previous EXOSAT observations (Wijnands & Reynolds 
2003). Therefore the quiescence timescale is currently 
unknown. The pulsar has a spin frequency of 314.4 
Hz and orbits in 4.3 hr around a ^ 0.1 Mq compan- 
ion (Markwardt & Swank 2003). We used all available 
X-ray data taken with the RXTE Proportional Counter 
Array in Event mode (time resolution 2~^^s) and Good 
Xenon mode data (time resolution 2^^°s) in the 2.5-15 
keV energy range. We closely follow the data reduc- 
tion procedure outlined in Watts et al. (2008) to which 
we refer for details. A total of 28 thermonuclear X-ray 
bursts are detected and removed in about 425 ks of data. 
Accretion-powered pulsations have a strong overtone at 
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twice the spin frequency (Watts et al. 2005). The time 
of arrivals (TOAs) of the pulsations were calculated by 
cross-correlating a ~ 500 s long folded profile with a sin- 
gle sinusoid whose frequency represents the NS spin fre- 
quency. The entire procedure is then repeated for the 
overtone. The effect of the Keplerian orbital motion is 
removed from the time of arrivals and a constant spin 
frequency is fitted to the TOAs. The X-ray flux is cal- 
culated as the total counts in each pulse profile, normal- 
ized to Crab units using the Crab spectral shape (van 
Straaten et al. 2003). 

2.2.1. The Phase-Flux Correlation 

XTE J1814 is the AMXP that shows the strongest lin- 
ear correlation between pulse phase and X-ray flux (Pa- 
truno, Wijnands & van der Klis 2009). This correlation 
has been interpreted as possibly being generated by a 
moving hot spot on the surface of the NS (Patruno, Wi- 
jnands & van der Klis 2009; Patruno et al. 2010). Such 
an interpretation is supported by recent results obtained 
by means of 3D MHD simulations of AMXPs (Romanova 
et al. 2003; Bachetti et al. 2010). Special care must be 
taken as flux induced pulse phase variations may resem- 
ble spin frequency variations. In Fig. 1 the effect of the 
phase-flux correlation is made evident by overplotting 
flipped pulse phase residuals (with respect to a constant 
spin frequency model) of the fundamental and the X-ray 
light curve. 
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Fig. 1. — Phase residuals (sign reversed) of the fundamental fre- 
quency (blue dots) and 2.5—16 kcV light curves ovorplottcd (red 
line) in arbitrary units. Phasc-fiux correlations can be seen on all 
timescales, with small deviations at the very end of the outburst. 
Similar results are obtained for the first overtone. 

It is obvious how variations of the pulse phase closely 
follow variations in the X-ray flux. In standard coher- 
ent pulsar timing the spin frequency derivatives are mea- 
sured by fitting a quadratic function to the pulse phases; 
by doing so for XTE J1814 the measured quadratic com- 
ponent will simply be a measure of the effect that the 
X-ray flux has on the pulse phases. 

2.2.2. Upper limits on the Spin Frequency Derivative 

To measure the spin frequency derivative we have first 
removed the effect of flux variations by using the lin- 
ear phase-flux correlation (for both harmonics). After 



"cleaning" the pulse phases, the residual rms variability 
is of the order of 0.02 cycles. We then fit a quadratic func- 
tion to the cleaned pulse phases and calculate the errors 
on the fitted spin parameters by means of ~ 10^ Monte 
Carlo simulations to account for the residual unmodclcd 
variability, as described in Hartman et al. (2008). 

The results indicate upper limits on the spin frequency 
derivative of the order of \v\ < 1.5 x 10~^^IIzs~^ at the 
95% confidence level. The upper limits are of the same 
order of magnitude as those reported by Hartman et al. 
(2008) for SAX J1808. These resuhs strongly contrast 
with the measurements obtained without cleaning the 
pulse phases (see e.g. Papitto et al. (2007) who report 

« -7 X lO-i-^Hzs-i). 

3. TORQUE ANALYSIS 

Standard accretion theory predicts a spin-up of the 
form (Bildsten et al. 1997): 

z> « 2.3 X 10-"C'/'Mrfo^''S8^'Af-f ^<^Hz/s, (1) 

where A/_io is the accretion rate in units of lO~^°M0/yr, 
i?8 the magnetic field in units of 10^ G, Mi. 4 the mass 
of the NS in units of IAMq, Riq its radius in units of 
10 km and ^ parametrises the uncertainties in evaluating 
the torque at the edge of the accretion disc and is thought 
to be in the range ^ « 0.3 — 1 (Psaltis & Chakrabarty 
1999). The average accretion rate for the outburst is 
~ 5 X lO-i°M0yr-i and 2 x 10-^°Mqyt-^ for SAX 
J1808 and XTE J1814, respectively, when considering 
the bolometric flux (using the data reported by Heinke et 
al. 2009; Wijnands & Reynolds 2003). The expression in 
(1) would thus give i> ~ 9 x lO^i-^Hzs-i for SAX J1808 
and !> ~ 4 lO^^Hzs^^ for XTE J1814, both larger than 
the observational upper limits. 

Our results suggest that during the outburst SAX 
J1808 and XTE J1814 spin up/down much less than can 
be measured and that the accretion torque is thus much 
weaker than the estimate in (1), or that there is an ad- 
ditional spin-down torque acting on the star, due for ex- 
ample to GW emission. 

3.1. Gravitational wave torques 

Gravitational wave emission was flrst suggested as the 
cause for the cutoff in the spin distribution of the LMXBs 
more than thirty years ago (Papaloizou & Pringle 1978). 
The main emission mechanisms that could be at work 
in these systems are crustal "mountains" (Bildsten 1998; 
Ushomirsky, Cutler & Bildsten 2000), magnetic deforma- 
tions (Cutler 2002; Melatos & Payne 2005) or unstable 
modes (Andersson 1998). All these processes can pro- 
duce a substantial quadrupole Q22 and thus a spin-down 
torque due to GW emission. 

3.1.1. Crustal mountains 

The crust of a NS is a thin layer (« 1 km thick) of 
matter in a crystalline phase, that can support shearing 
and the presence of a small asymmetry, or "mountain" . 
The crust consists of several layers of different nuclear 
composition and as accreted matter gets pushed further 
into the star it undergoes a series of nuclear reactions, 
including electron captures, neutron emission and pyc- 
nonuclear reactions (Sato 1979; Haensel & Zdunik 1990). 
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The heat deposited in the whole crust due to these reac- 
tions is Qh ~ 1-5 MeV per accreted baryon, but most of 
the heat is released by reactions occurring close to neu- 
tron drip (Haensel & Zdunik 1990). These reactions will 
heat the region by an amount (Ushomirsky & Rutledge 
2001): 

5T « 103C-V3o'Qn^A^22 K, (2) 

where Ck is the heat capacity in units oiks (the Boltz- 
mann constant) per baryon, pd is the pressure at which 
the reactions occur (in units of 10'^° erg/cm'^) and Qn 
the deposited heat per accreted baryon in MeV. AM is 
the accreted mass in units of 10^^ g, which for the sys- 
tems we consider is approximately the amount of mat- 
ter accreted after a few days of outburst. If the en- 
ergy deposition is (partly) asymmetric this would per- 
turb the equilibrium stellar structure and give rise to a 
mass quadrupole (Ushomirsky, Cutler & Bildsten 2000): 

fe-l-3xl0=x(^)(35|jv)'gc.„^ (3) 

where 5Tq is the asymmetric (quadrupolar) part of the 
temperature increase and Q is the threshold energy for 
the reaction. The quadrupole required for spin equilib- 
rium during an outburst is Q w 10'^^ g cm^ for both sys- 
tems (Ushomirsky, Cutler & Bildsten 2000). It is clear 
from equations (2,3) that even under the most optimistic 
assumptions it is very unlikely to build a "mountain" 
large enough to balance spin-up during accretion. 

3.1.2. Magnetic mountains 

It is well known that a magnetic star will not be spheri- 
cal and, if the rotation axis and the magnetic axis are not 
aligned, one could have a "magnetic mountain" leading 
to GW emission. However such deformations are unlikely 
to be large enough to balance the accretion torques in 
weakly magnetised systems such as the LMXBs (Haskell 
et al. 2008). Large internal fields could also cause a de- 
formation (Cutler 2002), but this would persist in qui- 
escence, leading to a rapid spin-down, of the order of 
the spin-up in (1), which is not observed in SAX J1808 
(Hartman ct al. 2009). 

Another possibility is that the magnetic field lines are 
stretched by the accreted material itself as it spreads 
on the star, giving rise to a large magnetically confined 
mountain. The results of Melatos & Payne (2005) sug- 
gest that the quadrupole built this way could balance the 
accretion torque only if the surface field is significantly 
stronger than the external dipole component. Further- 
more such a mountain would persist on an ohmic dissi- 
pation timescale r^^^ w 10^ yrs (Melatos & Payne 2005) 
and thus should also give rise to a strong additional spin- 
down in quiescence. 

3.1.3. Unstable modes 

A more promising scenario is that of an oscillation 
mode of the NS being driven unstable by GW emis- 
sion, the main candidate for this mechanism being the 
I = m ~ 2 r-mode(Andersson 1998). An r-mode is a 
toroidal mode of oscillation for which the restoring force 
is the Coriolis force. It can be driven unstable by the 
emission of GWs, as long as viscosity does not damp it 
on a faster timescale. This will only happen in a narrow 



window in frequency and temperature which depends on 
the microphysical details of the damping mechanisms (for 
a review see Andersson, Kokkotas & Ferrari (2001)). 

Of particular interest in this context is the situation 
in which the frequency at which the mode is driven un- 
stable (the 'critical' frequency) increases with tempera- 
ture in the range we are considering (T « 10^ — 10^ K). 
This could lead to the system being in equilibrium at 
the critical frequency and emitting GWs at a level that 
will balance the accretion torque (Andersson, Jones & 
Kokkotas 2002). This would be the scenario if the core 
contains hyperons (Nayyar & Owen 2006; Haskell & An- 
dersson 2010) or deconfined quarks (Andersson, Jones & 
Kokkotas 2002). If this is the case the dimensionless am- 
plitude of the mode a (as defined by Owen et al. 1998) 
will be approximately constant in time, i.e. a ~ 0, and if 
we assume that GW emission is balancing the accretion 
torque one has: 

« « 1.9 X 10-'M%lu-JJ'e^^Bl^\ (4) 

where t^4oo = J^/400 Hz and we have assumed a 1.4 Mq 
mass and 10 km radius, which we take as standard values 
from now on. 

After the outburst, when accretion ceases, the system 
will cool and spin-down following the critical frequency 
curve. Given that d « the mode amplitude would be 
approximately constant resulting in a spin-down rate ~ 
— 7 X 10^^^ Hz s~^ for both systems. This is considerably 
higher than the observed quiescent spin-down rate for 
SAX J1808 of j> « -5.5 x lO-^^ Hz s'^. 

Furthermore the shear from the mode will heat the star 
at a rate (Andersson, Kokkotas & Ferrari 2001): 

£: w 1.8 X lO^^a^rg'^/^i^logerg s-\ (5) 

where Tg = (T/10*A'). Let us first of all assume that the 
core cools due to the modified URCA neutrino emission 
process, which leads to an energy loss E sa 8.4 x lO'^^Tg 
erg/s (Page, Geppert & Weber 2006). Balancing this en- 
ergy loss with the viscous r-mode heating in (5) leads to 
an equilibrium temperature of T ~ 1.4 x 10* K for SAX 
J1808 and T « 1.5 x 10* K for XTE J1814. Let us now es- 
timate the temperature in the presence of the strongest 
possible neutrino emission processes. We therefore as- 
sume that nucleon direct URCA processes are possible, 
which may be the case in SAX J1808 (Yakovlev, Leven- 
fish & Haensel 2003), and that the core is not superfluid, 
with E^ ^ Ax lO^^Tl erg/s (Page, Geppert & Weber 
2006) (similar processes would be at work if hyperons or 
deconfined quarks are present). Note that in the pres- 
ence of superfluidity these processes would be strongly 
quenched (for a recent analysis of r-mode heating see Ho 
et al. (2011)). The temperatures we obtain for the stars 
are T « 1.5 x 10^ K for SAX J1808 and T « 1.6 x 10^ K 
for XTE J1814. 

How do these temperatures compare with observa- 
tional constraints? Both systems have upper limits on 
their temperature from quiescent observations (Heinke 
et al. 2009). One can map the surface temperature to 
the temperature at the top of the crust via the relation 
given by Potekhin, Chabrier & Yakovlev (1997) for an 
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accreted crust: 

[wkj =(lOi4cLs-0(^^-^10^j ' 

where g is the gravitational acceleration at the surface. 
By assuming that the star is isothermal, which is roughly 
correct if the thermal conductivity is high, as observa- 
tions of cooling X-ray transients suggest (Brown & Gum- 
ming 2009), one can then obtain a core temperature of 
T < 8.6 X 10^ K for SAX J1808 and T < 3.4 x 10^ for 
XTE J1814. 

It would thus appear that the presence of an unsta- 
ble r-mode is not consistent with the observed quiescent 
spin-down of SAX J1808, while it could be marginally 
consistent with the temperature of both sources only in 
the presence of direct URCA processes in the core. 

3.2. Accretion torques 

The interaction between the accretion disc and the 
magnetosphere is the natural candidate to explain not 
only the behaviour of our two systems, but also the cutoff 
in the spin-distribution of the millisecond pulsars. This 
mechanism was considered in detail by White & Zhang 
(1997), whose study seemed to indicated an unexpected 
link between the accretion rate and the magnetic field 
strength. This led to the suggestion that GW torques 
may be active, as in this case the dependence of the field 
strength on the accretion rate is weaker. 

Recent work has, however, shown that a more detailed 
modelling of magnetic torques in the disc can alleviate 
this problem (Andersson et al. 2005). It could, there- 
fore, be the case that both SAX J1808 and XTE J1814 
are close to the equilibrium spin rate set by such a pro- 
cess and that the torques are thus much weaker than 
the estimate in (1). In particular we will use the simple 
model of Andersson et al. (2005) to show that this is a 
more promising explanation than GW emission. 

We assume that a propeller phase starts at the reflaring 
stage for SAX J1808 (see a discussion of this possibility 
in Patruno et al. 2009). In XTE J1814 we assume that 
the onset of the propeller is at the point where the pulse 
phases deviate from a linear correlation (see Fig 1 and 
2). 

We calculate the average X-ray flux during an outburst 
and the flux at the onset of the propeller for XTE J1814 
(see Fig 2) and SAX J1808 and assume that the X-ray 
flux is a good tracer of the mass accretion rate. We 
assume that at the onset of the propeller the magneto- 
spheric radius (Bildstcn et al. 1997): 

r,„ = 35^M:',^'M-l/'Rll^'Bt/\m (7) 
is approximately equal to the co-rotation radius 

r, = 31A/i//i/4"oo/^km. (8) 

The results of (Andersson et al. 2005) indicate that the 
torque will vanish when the ratio between the propeller 
flux (fprop) and the average flux (favg) is such that 

{fpro-pl favgf'^ « 0.8. Wc find that {{fprop/favg?/^ « 

0.75 for XTE J1814 and 0.75-0.84 for SAX J1808 in the 
2002, 2005 and 2008 outbursts^ These values are very 

^ Only these outbursts had good enough observational coverage 
to allow the calculation. 



close to the theoretical value 0.8, which is encouraging 
given the large uncertainties in both the theoretical cal- 
culations and the observations. 

This does not mean that the accretion torques vanish 
during an outburst, but only that the spin variations can 
be smaller than those predicted by standard accretion 
theory. 
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Fig. 2.— X-ray lightcurve [2.5-16 koV] of XTE J1814. The 
dashed blue line indicates the possible onset of the propeller 
{rm = Tc) while the dashed black line is the average flux, which is 
close to the equilibrium condition r-m ~ 0.8 r'c. 



4. CONCLUSIONS 

We have presented a new analysis of the torques act- 
ing on XTE J1814 during the 2003 outburst and have 
interpreted the observed long term spin-down of SAX 
J1808 over a 10 yr baseline. We have found that there 
are no significant spin frequency derivatives that can be 
measured during the 2003 outburst of XTE J1814, with 
upper limits of I ;> I < 1.5x 10~^^ Hz s~^, of the same order 
of magnitude as the values determined for SAX J1808: 
\i>\ < 2.5 X 10-14 Hz s-i (Hartman et al. 2008, 2009). 
This implies that accretion torques during an outburst 
are smaller than expected from standard accretion the- 
ory for both AMXPs. 

We have analysed the two main mechanisms that could 
explain such a feature, namely the presence of GW 
torques and the possibility that the system may be close 
to spin equilibrium, as set by the disc/magnctosphere in- 
teraction. For SAX J1808 it is reasonably safe to say that 
GW torques can be excluded. For XTE J1814 a definite 
conclusion is harder to draw, as only one outburst has 
been observed and thus we do not have a measurement 
of the spin-down in quiescence. GW emission may still 
be marginally consistent with observations. 

A more promising possibility is that both sys- 
tems are close to spin equilibrium, as set by the 
disc/magnetosphere interaction. We have shown that the 
simple model of Andersson et al. (2005) predicts spin 
equilibrium close to the average accretion rate for the 
outbursts, thus leading to a reduced torque. Further 
work is necessary to analyse the problem with the use 
of refined disc models (e.g. D'Angelo & Spruit 2011), a 
task which is beyond the scope of this letter. 
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Note that the case of bright sources such as Sco X-1 
could be quite different, as high accretion rates would 
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